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ABSTRACT 

High-dynamic -range surface photometry in a companion paper makes possible accurate measurement of the 
stellar light deficits Ljef and mass deficits Mjef associated with the cores of elliptical galaxies. We show that Ldef 
correlates with the velocity dispersion a of the host galaxy bulge averaged outside the central region that may be 
affected by a supermassive black hole (BH). We confirm that Ldef correlates with BH mass M, . Also, the fractional 
light deficit Ljef/i'* correlates with M,/M^,, the ratio of BH mass to the galaxy stellar mass. All three correlations 
have scatter similar to or smaller than the scatter in the well known correlation between M, and a. The new 
correlations are remarkable in view of the dichotomy between ellipticals with cores and those with central extra 
light. Core light deficit correlates closely with M, and cr, but extra light does not. This supports the suggestion 
that extra light Es are made in wet mergers with starbursts whereas core Es are made in dry mergers. After dry 
mergers, cores are believed to be scoured by BH binaries that fling stars away as their orbits decay or by BHs that 
sink back to the center after recoiling from anisotropic gravitational radiation emitted when they merge. Direct 
evidence has been elusive. We interpret the new correlations as the "smoking gun" that connects cores with BHs. 
Together, the M, - a and M, - Ldet correlations give us two independent ways to estimate BH masses in core Es. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: formation — 
galaxies: nuclei — galaxies: photometry — galaxies: structure 



L INTRODUCTION 

Research on supermassive black holes (BHs) in galaxy nuclei 
centers on the growing evidence that BH growth and galaxy 
evolution are connected (Ho 2004). The most striking evidence 
is the discovery of correlations between BH mass M, and 
bulge (but not disk) luminosity (Kormendy 1993; Kormendy & 
Richstone 1995; Magorrian et al. 1998; Kormendy & Gebhardt 
2001) and between M, and the velocity dispersion cr of the 
bulge outside the region influenced by the BH (Ferrarese & 
Merritt 2000; Gebhardt et al. 2000b). These scaling relations 
are clues to formation physics. But they do not point to specific 
astrophysical processes. So the literature on their interpretation 
has been extraordinarily wide-ranging. 

This paper reports tight mutual correlations between M,, cr, 
and the light deficits Ldet that define cores in elliptical galaxies. 
These point more directly to formation processes. In particular, 
we find a correlation between the fundamental observables Ldet 
and cr that is as tight as the well known M, - a correlation. 

To set the stage, we recall that it is difficult to understand 
how cores can form in galaxies that are made by major mergers. 
We defin^ a core as the region near the center olF a bulge or 
elliptical interior to the radius where the outer, steep brightness 
profile shows a downward break to a shallow, central profile 
(Kormendy 1999; Kormendy et al. 2008). Caon et al. (1993) 
showed that the outer profiles I{r) are best modeled by Sersic 
(1968) functions, log/(r) oc r'/". Many authors confirm this. 
Kormendy et al. (2008) find that Sersic functions fit ellipticals 
with remarkable precision over large dynamic ranges. With 
respect to these outer Sersic fits, cores show central light 
deficits. Why? 

The problem is that galaxy mergers tend to preserve the 
highest progenitor densities. Lower-luminosity galaxies have 
higher central densities (Kormendy 1985, 1987; Lauer 1985; 



Faber et al. 1997). When ellipticals or bulges merge, this 
tends to destroy the above correlation (Faber et al. 1997). A 
possible solution is that cores may be scoured by the orbital 
decay of binary BHs that form in galaxy mergers Faber et al. 
1997; Milosavljevic & Merritt 2001; Milosavljevic et al. 2002; 
Merritt 2006). The orbit shrinks as the binary flings stars away. 
This decreases the surface brightness and excavates a core. The 
effect of a series of mergers should be cumulative; if the central 
mass deficit after one merger is /M,, then the mass deficit after 
major mergers should be ^ N fM, . The above papers predict 
that / ~ 0.5 to 2. Past observations of mass deficits Mdef were 
roughly consistent with this picture; Mdet oc M, and Nf ~ 1 - 5, 
consistent with formation by several successive dissipationless 
mergers (Milosavljevic et al. 2001, 2002; Ravindranath et al. 
2002; Graham 2004; Ferrarese et al. 2006; Merritt 2006). 

2. CORE MASS DEFICIT VERSUS BH MASS 

The key to this paper is high-accuracy surface photometry of 
all known ellipticals in the Virgo cluster presented in Kormendy 
et al. (2008, hereafter KFCB). Composite profiles are measured 
and assembled from 6-11 sources for each core elliptical. A 
variety of telescopes including the Hubble Space Telescope 
(HST) and wide-field, ground-based telescopes provide large 
dynamic range. Comparison of many data sources minimizes 
systematic errors. KFCB find that Sersic functions accurately 
fit 93 -99 % of the light of each galaxy. Central departures of 
the profiles from these fits can be interpreted with confidence. 
One result is an improved demonstration that profile shape 
participates in the well known dicotomy into two kinds of 
eflipticals (Bender 1988; Bender et al. 1989; Nieto, Bender, 
& Surma 1991; Kormendy & Bender 1996). HST photometry 
shows that ellipticals come in core and coreless varieties (Lauer 
etal. 1995, 2005, 2007; Gebhardt etal. 1996; Rest etal.2001; 
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Ravindranath et al. 2001). Faber et al. (1997) show that this 
distinction correlates with the rest of the dichotomy: core 
ellipticals generally have boxy-distorted isophotes and rotate 
slowly, while coreless ellipticals have disky-distorted isophotes 
and rotate rapidly. KFCB confirm results of Kormendy (1999) 
that coreless Es do not have featureless "power law profiles" 
near their centers but rather show distinct components - extra 
light above the inward extrapolation of their outer Sersic 
profiles. The extra light resembles the extra central components 
predicted in «-body simulations of galaxy mergers in which 
gas dissipation feeds starbursts (Mihos & Hernquist 1994; see 
Hopkins et al. 2008 for extensive simulations). KFCB conclude, 
in agreement with Faber et al. (1997), that "power law Es" (now 
more accurately called "extra light Es") formed in dissipative 
("wet") mergers with starbursts, whereas core Es formed in 
dissipationless ("dry") mergers with binary BH scouring. 

KFCB calculate with improved accuracy the mass excesses 
in extra light and the mass deficits (not) in cores. Figure 1 
shows that extra light correlates only weakly with M,. But 
mass deficits correlate closely with M,. Moreover, the mass 
deficit is a larger multiple of the BH mass than previously 
thought. The mean <logMdef/M,> = 1.02 ±0.07. These 
values are surprisingly large in comparison to the prediction 
(Merritt 2006) that Mdef/M, ~ 0.5 per major merger. However, 
with a more accurate treatment of the late stages of BH mergers, 
Merritt, Mikkola, & Szell (2007) find that Mdet /M, can be as 
large as 4 per merger. Second, an additional process has been 
proposed to make large-Mdet/M, cores (Merritt et al. 2004; 
Boylan-Kolchin, Ma, & Quataert 2004; Gualandris & Merritt 
2008). Coalescing binary BHs emit gravitational radiation 
anisotropically; they recoil at velocities comparable to galaxy 
escape velocities. If they do not escape, they decay back to the 
center by dynamical friction. In the process, they heat the core. 
Gualandris & Merritt (2008) estimate that they can excavate as 
much as Mjef/M, 5 in addition to the mass that was scoured 
by the pre-coalescence binary. So our observations present no 
problem for the idea that cores are made by a combination of the 
above BH scouring mechanisms acting over the course of one or 
more successive dry mergers. We do not estimate the number of 
mergers, because theoretical predictions of / are uncertain and 
because our measurements of Mjef depend on the assumption 
that unscoured profiles are exactly Sersic. 

The purpose of this Letter is to pursue the correlation in 
the bottom panel of Figure 1 in more detail. In particular, we 
"deconstruct" it to correlations between basic observables. 

First, an explanation of the measurements is necessary. The 
amount of extra or missing light is calculated by integrating 
the two-dimensional brightness distribution of the galaxy from 
the center to the inner limit rmin of the outer Sersic function's 
radial fit range. From this luminosity, we subtract the integral 
of the fitted Sersic function over the same radial range, keeping 
the ellipticity e of the Sersic function fixed at e(r,nin). Internal 
errors are estimated by substituting plausible extrapolations of 
the e(r) profile into the region at r < r^^^. Our measurements 
depend critically on the high accuracy and large dynamic range 
of the composite brightness profiles in KFCB; these constrain 
Sersic n values better than previous work. Our values of n 
are slightly larger than those derived previously (see KFCB), 
resulting in larger estimates of light deficits. Also, we calculate 
the difference between the integral of our outer Sersic fit and 
the integral of the galaxy profile, not the difference between a 
Sersic fit and a core-Sersic fit as in Ferrarese et al. (2006). 
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Fig. 1. — Stellar mass missing in cores {lower panel) or extra in coreless 
galaxies (upper panel) as a function of BH mass. Large and small symbols 
are for galaxies with and without dynamical BH detections. The sample is 
all known ellipticals in the Virgo cluster plus non-members denoted by open 
symbols. From Kormendy et al. (2008) with NGC 3379 added (logM. = 8.0). 

Missing light is converted to missing mass using mass-to-light 
ratios M/Ly oc L°^^ fitted to the SAURON sample of Cappellari 
et al. (2006) (zeropoint M/Ly = 6.07 at Mv = -21.6; estimated 
error in logM/Ly = 0.153). For galaxies with BH detections, 
we take M/Ly and M, from the dynamical modeling paper. 
Otherwise, M, is derived using the M, - a correlation as fitted 
by Tremaine et al. (2002). Velocity dispersions are mostly 
from Tremaine et al. (2002), Bender, SagHa, & Gerhard (1994), 
and Kormendy & Gebhardt (2001). The sample, total absolute 
magnitudes Myr, and light deficits My.dei are from KFCB with 
NGC 3379 added. Distances are fromMei et al. (2007) or Tonry 
et al. (2001). The data are listed in Table 1. 
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Core and BH Data 
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Fig. 2. — Fraction of the total V-band luminosity "missing" in cores versus 
the ratio of BH mass to the total stellar mass of the galaxy. The sample is as 
in Figure 1 and Table 1 . Large and small symbols denote galaxies with and 
without dynamical BH detections. The galaxies with BH detections are (top 
to bottom) M 87 (Macchetto et al. 1997), NGC 4261 (Ferrarese, Ford, & Jaffe 
1996), NGC 4374 (Bower et al. 1998), NGC 4649 (Gebhardt et al. 2003), and 
NGC 3379 (Gebhardt et al. 2000a). The line is Equation (1). 

3. CORRELATIONS BETWEEN LIGHT DEFICIT AND BH MASS 

Figure 2 shows the fraction of the total V-band light of the 
galaxy that is "missing" in the core versus the ratio of the BH 
mass to the stellar mass of the galaxy. As we would expect if 
Mdef/M, ~ constant (Figure 1), there is a good correlation and 
the slope is consistent with 1. 

The lines in Figures 2 - 4 are least-squares fits following 
Tremaine et al. (2002): (i) the dependent and independent 
variables x and y are treated symmetrically; (ii) we calculate 
regressions of y- j on x-x so errors in the fit parameter are 
essentially uncorrelated and minimized (Tremaine et al. 2002, 
p. 742); (iii) we calculate a fit using the estimated errors in both 
parameters and adding intrinsic scatter to bring the reduced 
to 1 ; (iv) we calculate a second fit by assigning equal errors to 
all X and different equal errors to all y and finding the values 
of those errors that result in a reduced x'^ of 1, thus providing 
partial protection against the danger that the assumed errors are 
not realistic, and (v) we adopt the mean of these fits as was done 
by Tremaine et al. (2002). In the above, e. g., x is the mean of 
the X values. In Figure 2, the adopted correlation is: 



/Lv 



0.01 



1.066: 



+0.236 
'-0.193 



M,/M 
0.001 



I.123±0.235 



(1) 



with RMS scatter = 0.23 in logLy^ef/Lv and 0.23 in logM./M. 
The latter is formally less than the error ±0.3 in logM, adopted 
by Tremaine et al. (2002). The fraction of the light that is 
missing in cores may be as good a predictor of M, as is a. 

The canonical mean mass fraction in BHs is M,/Mc:± 0.0013 
(Merritt & Ferrarese 2001; Kormendy & Gebhardt 2001). The 
5 BH detections in Figure 2 are consistent with this value. 
But if velocity dispersions reliably tell us BH masses, then 
there is significant scatter in M,/M toward lower values. Our 
observation that these lower M, values are correlated with 
lower Ldef values adds confidence to our conclusion that M,/M 



has significant cosmic scatter. The total range of logM,/M 
values found by Merritt & Ferrarese (2001) and by Kormendy 
& Gebhardt (2001) is more than ±1 around the mean. 

Figure 3 shows the absolute magnitude of the light that is 
missing in cores versus BH mass. It is analogous to previous 
derivations (see references in § 1) of the correlations of mass 
deficits with BH masses, but Mydef does not involve the 
uncertainty of mass-to-light ratios. 
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Fig. 3. — Absolute magnitude of the light that is missing in cores versus BH 
mass. The sample and symbols are as in Figure 2. The line is Equation (2). 



The correlation in Figure 3 is remarkably good. The fit is: 



Mv,dcf = (-18.05 ± 0.22) -(2.63 ± 0.42) log 



10^ Mp, 



(2) 



with RMS scatter = 0.52 in My^def and 0.20 in logM,. The latter 
is less than the error ±0.3 in logM, taken from Tremaine et al. 
(2002). Again, core light deficit may be as good a predictor of 
M, as is (T. Figure 3 spans the range of dynamically measured 
BH masses observed to date in ellipticals that have cores. 
For convenience, we can rewrite Equation (2) as: 



M. 
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108 L, 



V© 



0.95: 



40.18 
-0.13 



(3) 



Remarkably, the typical M/Ly ^ 8 almost cancels the factor-of- 
10 ratio between Ly^def/Lv and M,/M in Equation (1). Within 
the errors, M, ~ Ly^def, where both are measured in Solar units. 

We emphasize that the My.def measurements in Figure 3 are 
independent of M, except for the indirect dependence of M, 
on the light distribution through the dynamical modeling. It is 
interesting that the scatter in Figure 3, which mainly involves 
central galaxy quantities, is marginally smaller than the scatter 
in Figure 2, which also involves the total galaxy luminosity. 

4. LIGHT DEFICIT AND GALAXY VELOCITY DISPERSION 

The observation that M, correlates closely with both <t and 
My def hints that there may be a good correlation between My def 
and velocity dispersion. Figure 4 shows that this correlation is, 
in fact, excellent. Moreover, it is a correlation between pure 
observables, with no dynamical modeling involved. 
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Fig. 4. — Absolute magnitude of the light missing in cores versus 
galaxy velocity dispersion averaged inside a slit aperture whose length is 
approximately twice the half-light radius (Table 1). 

The correlation in Figure 4 is: 



Mvdef= (-16.71 ±0.18)-(13.73± 1.92) log, , 

, 250 km s 



108 L 



VQ 



= 4.10!!J 



63 



250 km s" 



5.49±0.77 



(4) 
(5) 



The RMS scatter is 0.55 in Mv.deU 0.22 in logLy.def, and 
0.0401 in logcr. Tremaine et al. (2002) estimate that the 
intrinsic dispersion in M, at constant a is 0.23 ± 0.05 dex. Our 
observed scatter in Ly def is essentially the same and includes 



measurement errors. For our sample, a predicts Lv_def slightly 
more accurately than it predicts M, . 

Figures 2 - 4 are remarkably good correlations, considering 
the difficulty of the def and M, measurements. The care that 
we took with the KFCB photometry helps. But there was never 
a guarantee that the intrinsic scatter in such correlations would 
be small. It is especially surprising because only core galaxies 
participate. Extra light galaxies contain BHs too (Figure 1), but 
extra light shows little correlation withM,. KFCB suggest that 
core Es are made in dry mergers whereas coreless Es are made 
in wet mergers in which central starbursts build the extra light. 
Then the amount of extra light is controlled mainly by the gas 
fraction in the progenitors and the properties of the starburst 
and only secondarily by the energy feedback from young stars, 
supernovae, and nuclear activity. We regard Figures 1-4 as 
further evidence for this explanation of the E-E dichotomy. 

In core galaxies, the excellent mutual correlations between 
light deficit, velocity dispersion, and BH mass - with observed 
scatter roughly consistent with measurement errors - imply 
that cores and BHs are closely related. We interpret them as 
a "smoking gun" in favor of the idea that cores are scoured 
by BHs. Core formation was surprisingly regular, given the 
probably heterogeneous merger histories of elliptical galaxies. 

Our "bottom line" conclusion is that, for the present galaxy 
sample, Ly def and a provide two independent and comparably 
accurate predictors of BH mass. 
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